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ABSTRACT 

We  investigate  the  relationship  of  the  H  n  region  luminosity  function  (H  n  LF)  to  the  H  n  region  size 
distribution  and  density-wave  triggering  in  grand-design  spiral  galaxies.  We  suggest  that  the  differential 
nebular  size  distribution  is  described  by  a  power  law  approximately  of  slope  —4,  with  flattening  at  radii  below 
~  1 30  pc.  This  contrasts  with  the  conventional  exponential  description,  but  it  is  physically  and  quantitatively 
consistent  with  the  typical  observed  value  of  -2  for  the  H  n  LF  slope.  To  study  H  n  LF  evolution,  we  have 
developed  an  interactive  code  that  computes  spatial  isochrones  for  the  evolving  loci  of  spiral  density  waves  in 
disk  galaxies.  This  allows  comparison  of  the  nebular  spatial  distribution  with  the  spatial  isochrones  for  sim¬ 
ple  rotation  curve  parameters.  Our  comparisons  for  four  grand-design  galaxies  suggest  that  the  corotation 
radius  rco  coincides  with  the  outer  ends  of  the  star-forming  arms.  This  value  for  rco  yields  the  best  spatial  cor¬ 
respondence  between  the  H  n  regions  and  the  isochrones  and  also  appears  to  yield  a  coincidence  between  the 
inner  Lindblad  resonance  with  the  radial  onset  of  star  formation  in  the  arms.  Thus,  we  suggest  that  isochrones 
offer  a  new,  simple,  and  effective  technique  for  determining  rco  and  thus  the  spiral  pattern  speed.  However, 
application  of  the  isochrones  also  demonstrates  that  the  evolution  of  the  nebular  population  is  difficult  to 
spatially  isolate  in  these  galaxies. 
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1.  INTRODUCTION 

The  H  ii  region  luminosity  function  (H  n  LF)  is  an 
important  diagnostic  of  star  formation  properties  in  gal¬ 
axies.  It  is  a  fairly  straightforward  parameter  to  determine 
for  nearby  galaxies  (e.g.,  Kennicutt,  Edgar,  &  Hodge  1989; 
Banfi  et  al.  1993;  Rozas,  Beckman,  &  Knapen  1996a; 
Thilker  et  al.  2002)  and  is  usually  presented  in  the 
differential  form  as  a  power  law: 

N(L)  dL  =  AL~a  dL  ,  (1) 

where  N(L )  dL  is  the  number  of  nebulae  with  Ha  luminosi¬ 
ties  in  the  range  L  to  L  +  dL.  Occasionally,  the  cumulative 
H  ii  LF  is  given,  but  we  note  that  the  differential  form  is  a 
more  powerful  diagnostic  of  the  underlying  properties 
(Cioffi&  Shull  1991). 

In  an  earlier  paper  (Oey  &  Clarke  1998,  hereafter  Paper 
I),  we  modeled  the  behavior  of  the  H  n  LF  and  identified 
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ways  in  which  it  can  be  used  to  examine  the  current  proper¬ 
ties  of  massive  star  formation  and  recent  (<10  Myr)  star 
formation  history.  We  assumed  a  power-law  distribution  in 
A*,  the  number  of  ionizing  stars  per  object, 

A(A* )  dN*  =  BN dN*  .  (2) 

We  then  performed  Monte  Carlo  simulations,  in  which  we 
draw  A*  from  this  distribution  and  also  draw  the  individual 
stellar  masses  from  a  Salpeter  (1955)  initial  mass  function 
(IMF).  These  models  showed  that  reported  variations  of  the 
measured  slope  a  (eq.  [1])  in  different  environments  do  not 
require  varying  /?,  the  slope  of  the  A*  distribution  (eq.  [2]) 
but,  rather,  can  be  explained  by  variations  in  the  star  forma¬ 
tion  history  and  the  upper  cutoff  in  the  A*  distribution.  An 
important  feature  is  that  stochastic  effects  for  objects  with 
small  A*  cause  a  flattening  in  the  H  n  LF  below  a  certain 
luminosity  Lsat.  Sparse  clusters  subject  to  this  effect  are 
termed 44  unsaturated  ”  with  respect  to  the  IMF,  and  clusters 
with  higher  A*  we  call 44  saturated.”  We  identified  Lsat  with 
the  Ha  luminosity  associated  with  the  most  massive  star  in 
the  IMF,  although  we  refer  readers  to  Paper  I  for  a  complete 
discussion. 

A  basic  test  on  the  nature  of  the  H  n  LF  is  its  relation  to 
the  size  distribution  of  the  H  n  regions.  Although  the  H  n 
LF  has  been  empirically  described  by  a  power  law,  the 
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Fig.  1 . — Model  H  n  LFs:  (a)  zero-age  and  ( b )  evolved  single-burst,  from  Oey  &  Clarke  (1998,  Figs.  2 a  and  2c).  The  dotted  line  marks  Lsat. 


corresponding  size  distribution  has  been  described  by  an 
exponential  law  (van  den  Bergh  1981;  Hodge  1983,  1987). 
As  discussed  above,  Paper  I  found  that  the  slope  breaks  and 
other  features  in  the  H  n  LF  are  nevertheless  best  explained 
by  a  simple  power  law  in  A*  (eq.  [2]  of  Paper  I).  Therefore, 
we  should  naively  expect  the  size  distribution  to  also  follow 
a  power-law  distribution.  In  §  2  below,  we  examine  whether 
observations  of  the  H  n  region  size  distribution  are 
consistent  with  a  simple  power-law  relationship. 

Another  test  on  the  nature  of  the  H  n  LF  addresses  its 
evolutionary  properties.  The  models  constructed  in  Paper  I 
included  the  evolution  of  the  H  n  LFs  resulting  from  a  single 
burst  creation  scenario.  In  Figure  1,  we  show  representative 
zero-age  and  7  Myr  models  (Figs.  2 a  and  2c  of  Paper  I).  The 
zero-age  model  shows  the  flattening  below  Lsat.  As  time  goes 
by,  the  entire  H  n  LF  shifts  to  lower  L.  However,  when  con¬ 
sidering  a  fixed  range  of  L ,  this  can  result  in  an  apparent 
elimination  of  the  slope  break  in  strongly  evolved  popula¬ 
tions,  due  to  the  differential  fading  effect.  If  this  model  turns 
out  to  be  correct,  it  would  then  be  possible  to  distinguish, 
for  example,  between  an  older  nebular  population,  in  which 
high-luminosity  objects  were  created  in  a  single  burst,  and 
the  continuous  creation  of  low-luminosity  H  n  regions. 

Paper  I  applied  these  results  to  observations  of  arm  and 
interarm  regions  of  grand-design  spiral  galaxies,  suggesting 
that  arm  populations  represent  the  current  burst  of  star 
formation,  while  interarm  objects  represent  an  aged  popula¬ 
tion  remaining  in  the  wake  of  the  spiral  density  waves. 
Observations  in  the  literature  sometimes  suggest  that  inter¬ 
arm  populations  show  a  slightly  steeper  H  n  LF  than  the 
arm  populations  (e.g.,  Kennicutt  et  al.  1989;  Rand  1992). 
This  phenomenon  can  be  explained  by  the  two-slope 
structure  in  the  zero-age  models,  which  will  cause  a  slope 
measurement  to  be  flatter  when  fitted  over  L  below  Lsat  (see 
Fig.  1).  In  contrast,  the  evolved  populations  show  the  intrin¬ 
sic  slope  of  the  high-L  regime.  Figure  5  of  Paper  I  shows  that 
this  interpretation  of  evolution  between  the  arm  and 
interarm  H  n  regions  is  consistent,  in  all  cases,  with 
available  observations  of  six  grand-design  spiral  galaxies. 

Observations  of  the  H  n  LF  in  the  literature  to  date  have 
been  subjectively  and  coarsely  binned  into  the  arm  and 
interarm  categories.  A  more  powerful  test  of  this  evolution¬ 
ary  scenario  would  be  to  rebin  the  H  n  regions  more  quanti¬ 
tatively  into  regions  corresponding  to  fixed  age  ranges 
behind  the  leading  edge  of  the  spiral  arms,  as  inferred  from 
the  spiral  arm  pattern  speeds.  We  might  then  be  able  to 


examine  the  competing  contributions  in  interarm  regions  of 
aging  nebular  populations  versus  bona  fide  interarm  star 
formation.  These  are  our  motivations  for  carrying  out  the 
studies  described  in  §  3.  As  an  important  by-product  of  this 
study,  we  developed  a  technique  that  compares  nebular 
positions  with  spatial  isochrones  for  the  spiral  density 
waves.  This  offers  a  new  method  for  locating  the  corotation 
resonances  in  spiral  galaxies  and  is  described  in  §  4. 

2.  THE  FORM  OF  THE  H  n  REGION 
SIZE  DISTRIBUTION 

As  mentioned  in  §  1,  the  conventional  parameterization 
of  the  H  ii  region  size  distribution  as  an  exponential  (e.g., 
van  den  Bergh  1981;  Hodge  1983, 1987)  is  intuitively  at  odds 
with  the  power-law  form  of  the  H  n  LF.  Since  the  behavior 
of  the  H  ii  LFs  appears  to  be  fairly  well  understood  (e.g., 
Paper  I),  it  is  important  to  reexamine  the  size  distribution 
within  the  same  framework,  since  they  are  physically 
related.  In  particular,  we  might  also  expect  the  size 
distribution  to  follow  a  power  law  of  the  form 

Nr(R)  dR  oc  R~b  dR  ,  (3) 

where  R  is  the  nebular  radius.  This  has  indeed  been 
suggested  for  some  galaxies  (e.g.,  Pleuss,  Heller,  &  Fricke 
2000). 

Since  the  nebular  Ha  luminosity  L  is  integrated  from  the 
nebular  volume  emission,  we  should  expect  L  oc  R3  and 
dL/dRoc  R2.  Therefore,  the  differential  size  distribution 
should  be 

Nr(R)  dR  =  N(L)  ^  dR  oc  R2~2a  dR  .  (4) 

dR 

Hence  we  find  that  the  slopes  a  and  b  of  the  H  n  LF  and  size 
distribution,  respectively,  are  related  as  b  =  2  —  3a.  For  the 
usual  value  a  =  2  found  for  the  high-luminosity,  44  satu¬ 
rated”  H  ii  regions,  we  therefore  predict  a  corresponding 
size  distribution  having  b  =  4. 

As  mentioned  earlier,  a  turnover  in  the  H  ii  LF  caused  by 
stochastic  effects  for  44  unsaturated  ”  objects  occurs  at  a  lumi¬ 
nosity  Lsat.  Since  Lsat  is  associated  with  the  most  massive 
stars  in  the  IMF,  their  Stromgren  radius  Rsat  should  be  a  cor¬ 
responding  inflection  point  in  the  size  distribution.  This  flat¬ 
tening  at  small  R  can  mimic  an  exponential  distribution.  For 
the  most  massive  stars,  stellar  atmosphere  models  suggest 
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Fig.  2. — Differential  H  n  region  diameter  distributions  for  galaxies  in  Table  1,  fitted  with  power  laws  for  logD/pc  >2.3.  The  fitted  values  of  —b  and  the 
standard  deviation  are  shown  [note  that  the  fitted  slope  for  log  7V(log  D)  vs.  log  D  is  1  -  b\.  Empty  bins  are  shown  as  log  7V(log  D)  =  - 1 . 


logLsat/ergs  s-1  ~  38.0-38.5  (Schaerer  &  de  Koter  1997; 
Vacca,  Garmany,  &  Shull  1996;  Panagia  1973),  implying 
Rsat  ~  130  pc  for  a  density  of  1  cm-3.  This  estimate  for  Rsat 
should  be  taken  as  an  upper  limit,  since  the  density  is  likely  to 
be  higher  in  the  actual,  central  star-forming  area. 

We  can  now  compare  observations  of  the  H  n  region  size 
distribution  to  our  predictions  for  slope  b  =  4,  with  a  flat¬ 
tening  in  R  below  logRsat/pc  ~  2.1.  We  refitted  the  nebular 
diameter  distributions  with  power  laws  for  a  number  of  gal¬ 
axies  in  the  literature.  Since  size  determinations  can  be  sensi¬ 
tive  to  different  methods,  and  thus  to  ensure  relatively 
consistent  comparisons,  we  used  data  from  only  two  groups 
of  investigators:  the  IAC-Hertfordshire  collaboration 
(Knapen  et  al.  1993;  Rozas  et  al.  1996b;  Knapen  1998)  and 
Youngblood  &  Hunter  (1999).  Hodge’s  group  (e.g.,  Hodge 
1983,  1987;  Hodge  et  al.  1994,  1999)  and  Petit’s  group  (e.g., 
Petit,  Sivan,  &  Karachentsev  1988;  Courtes  et  al.  1993;  Petit 
1998)  have  contributed  many  studies  of  H  n  region  size  dis¬ 
tributions;  however,  their  target  galaxies  generally  do  not 
offer  enough  objects  to  adequately  sample  the  nebular  size 
distribution  at  log  D/ pc  >  2.4,  where  D  is  the  nebular  diam¬ 
eter.  Likewise,  we  include  only  the  three  galaxies  from  the 
Youngblood  &  Hunter  (1999)  sample  that  provide  adequate 
sampling  in  this  regime. 


Figure  2  shows  the  size  distributions  for  the  nine  galaxies 
listed  in  Table  1.  We  fitted  power-law  slopes  for  log  N(D)/ 
pc  >  2.3,  with  the  bins  weighted  by  y/N.  The  IAC- 
Hertfordshire  galaxies  are  grand-design  spirals  and  the 
Youngblood  &  Hunter  galaxies  are  Magellanic  irregulars. 

TABLE  1 

References  for  Galaxy  Diameter  Distributions 


Galaxy  Typea  Reference 


Ml  00 .  Sc(s)  I  1 

NGC157 .  Sc(s)  I— II  2 

NGC1156 .  Sm  IV  3 

NGC2366 .  SBmIV-V  3 

NGC3631 .  Sbc(s)  II  2 

NGC4214 .  SBmlll  3 

NGC6814 .  Sbc(rs)  I— II  4 

NGC6764 .  SBbc(s)b  2 

NGC6951 .  Sb/SBb(rs)  1.3  2 


a  Type  from  Sandage  &  Tamann  1987. 
b  Type  from  Rozas  et  al.  1996b. 

References. — (1)  Knapen  1998;  (2)  Rozas  et  al. 
1996b;  (3)  Youngblood  &  Hunter  1999;  (4)  Knapen 
et  al.  1993. 
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Both  sets  of  galaxies  have  H  n  LFs  extending  to  high  lumi¬ 
nosities,  consistent  with  no  upper  cutoff  in  luminosity.  Our 
fits  to  the  size  distributions  therefore  should  not  be  influenced 
by  possible  truncation  at  the  upper  end  of  the  H  n  LF.  It 
appears  that  the  H  n  region  diameter  distributions  shown  in 
Figure  2  are  consistent  with  power  laws  in  the  given  regime. 
Objects  with  smaller  diameters  show  a  flattening  in  the  size 
distribution,  qualitatively  consistent  with  the  flattening  in  the 
H  ii  LF  for  statistically  unsaturated  objects. 

The  small  dynamic  range  and  small  number  statistics 
cause  some  ambiguity  and  large  uncertainties  in  the  fitted 
slope  —b,  as  shown  in  Figure  2.  We  also  caution  that  we  can¬ 
not  rule  out  the  existence  of  systematic  effects.  Nevertheless, 
the  average  slope,  weighted  inversely  by  the  uncertainties,  is 
(b)  =  4.12  with  a  standard  deviation  of  0.90.  This  prelimi¬ 
nary  result  agrees  with  the  value  predicted  above  of  b  =  4. 
Figure  2  also  shows  that  the  turnover  in  the  diameter  distri¬ 
bution  is  also  consistent  with  the  predicted  value  of 
logD/pc  ~  2.4.  We  therefore  suggest  that  the  H  n  region  size 
distribution  is  intrinsically  described  by  a  power-law  function 
in  the  regime  beyond  the  turnover,  consistent  with  its  physi¬ 
cal  relation  to  the  H  n  LF.  While  consistent  functional 
forms  are  expected  for  the  size  distribution  and  H  ii  LF,  fur¬ 
ther  confirmation  of  the  empirical  power  law  for  the  former 
is  desirable,  with  attention  to  systematic  uncertainties  in  size 
determinations. 

3.  TESTING  FOR  AN  EVOLVING  H  n  LF 

Based  on  Paper  I,  we  approximate  that  the  H  n  regions 
found  along  the  leading  edge  of  spiral  arms  represent  the 
youngest,  zero-age  population  of  objects.  H  n  regions 
behind  the  leading  edge  therefore  represent  a  somewhat 
older  population  remaining  in  the  wake  of  the  spiral  density 
wave.  If  the  galactic  rotation  curve  and  spiral  arm  pattern 
speed  are  known,  it  is  in  principle  possible  to  estimate  the 
ages  of  objects  based  on  their  location  with  respect  to  the 
leading  edges  of  the  spiral  arms.  This  would  also  assume 
that  the  velocities  of  the  H  n  regions  do  not  deviate  signifi¬ 
cantly  from  the  galactic  rotation  curve  during  the  passage  of 
the  density  wave. 

We  have  developed  a  code  that  computes  spatial  iso¬ 
chrones  from  the  rotation  curve,  given  the  current  locus  of 
the  spiral  arms.  We  defer  discussion  of  the  code  and  spiral 
arm  kinematics  until  §  4  below,  and  here  use  the  spatial  iso¬ 
chrones  to  search  for  evolutionary  features  as  predicted  in 
Paper  I.  We  examined  four  of  the  grand-design  spiral  gal¬ 
axies  studied  in  Paper  I,  while  the  two  remaining  galaxies 
included  in  that  paper,  NGC  3631  and  NGC  6814,  did  not 
have  a  simple  two-arm  structure  suitable  for  this  analysis. 
Figure  4  below  shows  Ha  images  of  the  galaxies:  M51 
(Rand  1992),  Ml 00  (Knapen  et  al.  1993),  NGC  157,  and 
NGC  6951  (Rozas  et  al.  1996a).  Five  sets  of  spatial  iso¬ 
chrones  are  calculated  for  each  of  the  two  dominant  arms  in 
these  galaxies,  calculated  at  intervals  of  12  Myr.  The  zero- 
age  isochrones  correspond  to  an  actual  spatial  fit  to  the  arm 
positions.  The  adopted  corotation  radius  is  shown  by  the 
dashed  black  line,  and  the  inner  Lindblad  resonance  (ILR) 
and  outer  Lindblad  resonance  (OLR)  are  shown  by  the  solid 
and  dotted  black  lines,  respectively,  assuming  m  =  2  spiral 
arms  for  these  galaxies. 

Unfortunately,  it  appears  that  the  isochrone  spatial  sepa¬ 
ration  that  is  sensitive  to  the  nebular  evolutionary  timescale 
is  on  the  order  of  the  spatial  scatter  of  the  objects  along  the 


spiral  density  wave.  Figures  4 a-Ac  show  that  the  best- 
defined  arms  are  enclosed  within  the  —12  to  12  Myr  iso¬ 
chrones,  but  show  a  spatial  scatter  that  precludes  further 
distinction.  NGC  6951  is  an  exception  (Fig.  Ad),  showing 
larger  separation  between  the  isochrones  that  also 
corresponds  to  the  nebular  spatial  distribution. 

Figure  3  shows  the  H  n  LFs  for  subgroups  of  H  n  regions 
binned  within  isochrones  for  — 12  to  12  Myr  (top),  0-24  Myr 
(middle),  and  12-36  Myr  (bottom).  As  discussed  above  from 
Paper  I,  the  H  n  LF  naturally  shows  a  break  in  slope  below 
the  high-luminosity  power-law  tail,  which  is  caused  primar¬ 
ily  by  the  transition  to  unsaturated  stellar  clusters  domi¬ 
nated  by  small  number  statistics.  For  older  nebular 
populations,  this  slope  transition  at  Lsat  is  expected  to  shift 
toward  lower  luminosities,  while  the  high-luminosity 
power-law  slope  remains  the  same  (Fig.  1;  see  also  Paper  I). 
The  vertical  dotted  line  in  these  figures  shows  the  luminosity 
bin  in  which  we  estimate  Lsat  to  occur  for  the  youngest  bin¬ 
ning;  this  value  is  plotted  for  reference  and  remains  constant 
between  the  three  H  n  LF  subsets  for  each  galaxy. 

It  is  apparent  in  Figure  3  that  the  location  of  Lsat  is 
essentially  indistinguishable  between  the  three  subsets. 
The  transition  point  generally  appears  to  be  more  ambig¬ 
uous  in  the  oldest  binnings  (bottom)  as  compared  to  the 
youngest  (top),  but  more  galaxies  must  be  examined  to 
determine  whether  that  pattern  is  real.  If  so,  then  it 
would  be  consistent  with  a  greater  contribution  from 
evolved  populations  (Paper  I).  The  ambiguities  again 
reflect  our  discussion  above,  that  the  isochrones  for  these 
galaxies  are  insensitive  to  aging  effects,  since  the  spatial 
scatter  is  on  the  order  of  the  separation  for  isochrones 
over  the  nebular  evolution  timescale. 

A  number  of  other  factors  probably  also  contribute  to  the 
insensitivity  of  the  isochrones  to  aging  effects.  First,  the 
models  in  Paper  I  (e.g.,  Fig.  1)  assumed  the  instantaneous 
creation  of  all  objects  in  the  H  n  LF,  whereas  triggering  by 
the  spiral  density  wave  most  likely  has  a  finite  timescale. 
Since  most  H  n  regions  recombine  on  timescales  of  <10 
Myr,  we  would  not  expect  any  significant  population  to  be 
visible  at  24  Myr.  For  the  most  part,  this  is  consistent  with 
the  fairly  negligible  Ha  emission  at  the  24  Myr  isochrone 
for  most  of  the  arms  in  M51,  M100,  and  NGC  157.  NGC 
6951  shows  spatially  more  extensive  star  formation,  with 
H  ii  regions  also  seen  along  the  24  Myr  isochrone;  this  there¬ 
fore  implies  an  extended  timescale  for  star  formation  at  any 
given  location.  An  extended  star  formation  timescale  is  also 
equivalent  to  a  spatial  scatter  for  objects  of  any  given  age,  as 
mentioned  above.  For  the  H  n  LFs,  this  will  blur  the  distinc¬ 
tion  between  the  youngest  subset  and  the  next  to  youngest. 
It  may  also  be  that  the  luminosity  fading  of  the  H  n  regions 
is  even  faster  than  our  extreme  model  in  Paper  I  (see  Paper 
I,  Fig.  2f).  Another  factor  that  affects  the  quantitative  com¬ 
parison  to  the  H  ii  LF  evolution  is  the  uncertainty  in  identi¬ 
fying  and  fitting  the  locus  of  the  zero-age  density  wave. 
Finally,  our  simple  analytic  description  for  the  rotation 
curve  and  pattern  speed  are  idealized  forms  that  do  not 
describe  the  details  of  the  gas  kinematics  (see  §  4  below), 
which  also  affect  the  location  of  the  spatial  isochrones. 

Thus,  our  test  of  the  evolutionary  features  in  the  H  n  LF 
for  these  galaxies  is  inconclusive,  though  an  upper  limit  on 
the  evolutionary  effects  is  indicated.  For  most  of  the 
galaxies,  the  spatial  isochrones  do  not  offer  the  necessary 
separation  to  isolate  populations  of  H  n  regions  triggered  in 
the  arm  density  waves. 
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Fig.  3. — H  n  LFs  for  objects  binned  between  isochrones  for  ( a )  M51,  ( b )  M100,  (c)  NGC  157,  and  ( d )  NGC  6951.  Top ,  Bins  1  and  2  (—12  to  12  Myr); 
middle ,  bins  2  and  3  (0-24  Myr);  and  bottom,  bins  3  and  4  (12-36  Myr). 
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4.  LOCATING  COROTATION  WITH 
SPATIAL  ISOCHRONES 

We  now  describe  our  analysis  of  spatial  isochrones,  which 
offer  some  interesting  insights  on  spiral  arm  kinematics. 

While  measurements  of  rotation  curves  are  fairly  straight¬ 
forward,  historically  it  has  been  much  more  difficult  to  esti¬ 
mate  the  spiral  density-wave  pattern  speed  or,  equivalently, 
to  determine  the  location  of  the  corotation  radius.  Several 
methods  for  determining  the  corotation  radius  or  pattern 
speed  have  been  proposed  in  the  past.  The  earliest  proposal 
was  based  on  Lin’s  (1970)  conjecture  that  organized,  two¬ 
armed  spiral  patterns  exist  only  within  the  corotation  circle 
and  that  arms  become  more  fragmented  outside.  Shu, 
Stachnik,  &  Yost  (1971)  were  among  the  first  to  apply  this 
conjecture  to  physical  galaxies,  placing  the  corotation  where 
the  distribution  of  H  n  regions  is  seen  to  end.  Tremaine  & 
Weinberg  (1984)  developed  an  approach  for  pattern  speed 
determination  based  on  measurement  of  surface  brightness 
and  velocity  distributions  of  luminosity  density  tracers, 
assuming  continuity  in  the  tracer.  Canzian’s  (1993)  method 
for  determining  corotation  relies  on  the  qualitative  differ¬ 
ence  in  the  residual  velocity  field  inside  versus  outside  of 
corotation.  Elmegreen,  Elmegreen,  &  Seiden  (1989)  identi¬ 
fied  resonances  through  breaks  in  spiral  arms,  and  Cepa  & 
Beckman  (1990)  utilized  the  radial  star  formation  efficiency. 
Elmegreen,  Wilcots,  &  Pisano  (1998)  determined  corotation 
from  the  transition  of  the  inward-to-outward  streaming  gas 
velocity.  Finally,  Clemens  &  Alexander  (2001)  used  high- 
resolution  spectroscopic  data  to  find  signatures  of  the  pre- 
and  postshock  gas  across  the  spiral  arms;  this  information 
determines  the  shock  velocity  and,  hence,  the  pattern  speed. 
In  what  follows,  we  explore  a  new  approach  that  uses  the 
nebular  spatial  distribution  to  constrain  the  corotation 
radius. 

For  a  galaxy  with  a  fiat  rotation  curve,  the  corotation 
radius  rco  is  related  to  the  angular  pattern  speed  Vtp  as 


where  vc[rc  is  the  circular  velocity.  Assuming  that  star  forma¬ 
tion  is  triggered  by  spiral  density  waves,  the  locus  of  the 
youngest  H  n  regions  essentially  indicates  the  location  of 
the  density  wave.  As  the  nebulae  age,  their  positions  will  be 
determined  by  both  the  spiral  arm  pattern  speed  and  the  gal¬ 
axy’s  rotation  curve.  It  is  therefore  possible  to  determine 
spatial  isochrones  for  a  given  galaxy,  which  represent  the 
positions  of  objects  that  all  formed  simultaneously  within  a 
density  wave. 


Our  code  computes  spatial  isochrones,  given  an  initial 
locus  for  the  spiral  arm  pattern.  The  H  n  region  positions 
are  deprojected  and  converted  to  polar  coordinates,  and  we 
fit  the  spiral  pattern  at  the  arm’s  leading  edge,  using  a  poly¬ 
nomial  function.  We  then  compute  spatial  isochrones  from 
the  locus  of  the  spiral  arm  pattern,  rotation  curve,  and  puta¬ 
tive  corotation  radius.  We  assume  a  flat  rotation  curve, 
except  for  the  innermost  regions  r  <  rc irc,  where  we  assume 
v  oc  r  (Table  2). 

Our  code  includes  a  widget  to  easily  adjust  the  assumed 
corotation  radius.  Varying  rco  by  hand,  we  found  a  value  of 
rco  that  determines  a  set  of  isochrones  having  a  noticeable 
correspondence  to  the  spatial  configuration  of  the  H  n 
regions.  These  values  of  rco  turn  out  to  correspond  to  the 
ends  of  the  spiral  arms  in  the  star-forming  disk.  As  we 
discussed  above,  other  investigators  (e.g.,  Shu  et  al.  1971; 
Rohlfs  1977;  Kenney  et  al.  1992;  Zhang,  Wright,  Alexander 
1993)  have  in  the  past  identified  rco  with  the  ends  of  the 
spiral  arms.  For  three  of  our  four  galaxies,  Figure  4  shows 
that  these  values  of  rco  also  yield  ILRs  that  coincide  well 
with  the  radial  onset  of  significant  star  formation  in  the 
arms.  Figure  4b  also  suggests  a  correspondence  between  the 
OLR  and  the  edge  of  the  star-forming  disk  in  Ml  00. 

For  all  four  galaxies,  the  regions  near  the  ILRs  display 
the  characteristic  44  twin-peaks  ”  morphology  (Kenney  et  al. 
1992)  or  44  twin-peaks-plus-inner-bar  ”  (Zhang  et  al.  1993), 
together  with  an  annular  zone  between  the  twin  peaks  and 
the  inner  bar/nucleus  that  is  swept  clear  of  star  formation. 
This  combination  of  central  morphologies  is  common 
among  mid-  to  early-type  galaxies  and  is  related  to  the  orbi¬ 
tal  behavior  of  gas  near  the  double  ILRs  that  are  present  in 
these  galaxies  with  high  central  potential  wells  (Zhang  et  al. 
1993  and  references  therein),  namely,  the  inner  ILR  (IILR) 
and  the  outer  ILR  (OILR).  This  morphology  allows  the 
effective  central  channeling  of  gas  and  fueling  of  nuclear 
activity.  Note  that  our  simplified  assumed  rotation  curves 
yield  only  one  ILR,  which  presumably  is  located  between 
the  IILR  and  OILR  of  the  actual  rotation  curves. 

For  M51  and  Ml 00,  Figure  4  shows  that  setting  rco  to  the 
tips  of  the  spiral  arms  causes  the  isochrones  to  bunch  or 
converge,  but  not  cross  over,  at  azimuthal  locations  near 
large  clumps  of  nebulae.  These  clumps  are  not  confined 
within  the  isochrones  and  extend  to  radii  substantially 
beyond  them;  however,  their  azimuthal  coincidence  with 
the  isochrone  convergences  strongly  suggests  that  these  loci 
are  associated  with  dynamical  triggering  of  the  star  forma¬ 
tion.  These  patterns  are  most  likely  related  to  noncircular 
motions  and  asymmetry  in  the  mass  distributions.  Indeed, 
both  M51  (see  Garcia-Burillo  et  al.  1993b)  and  Ml 00 


TABLE  2 

Galaxy  and  Resonance  Parameters 


Distance  vcirc  rcirc  rco  ILR  OLR 


Galaxy  (Mpc)  (kms  ])  (kpc)  (kpc)  (kpc)  (kpc)  Referencea 


M51  .  9.6  220  0.93  14  4.1  24.9  1 

M100 .  16.1  210  1.56  12  3.5  20.5  2 

NGC  157  .  22.2  170  0.54  6  1.8  10.2  3 

NGC6951 .  19.0  235  0.46  10  2.9  17.1  4 


a  References  for  rotation  curve  parameters:  rcirc  is  the  radius  beyond  which  the  circular  velocity  is 
taken  to  be  constant,  at  value  Ecirc. 

References. — (1)  Rand  1993;  (2)  Sempere  et  al.  1995;  (3)  Sempere  &  Rozas  1997,  Ryder  et  al.  1998; 
(4)  Rozas  et  al.  2002. 
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Fig.  4  a 

Fig.  4. — Ha  images  of  the  sample  galaxies  overlaid  with  spatial  isochrones  for  —12,  0,  12,  24,  and  36  Myr.  Also  shown  are  the  adopted  corotation  radius 
(dashed  line),  the  ILR  (solid  line),  and  the  OLR  (dotted  line),  (a)  M51,  with  Ho  data  from  Rand  (1992).  (b)  M100,  with  Ho  data  from  Knapen  et  al.  (1993).  (c) 
NGC  157,  with  Ho  data  from  Rozas  et  al.  (1996a).  (d)  NGC  6951,  with  Ho  data  from  Rozas  et  al.  (1996a)  [See  the  electronic  edition  of  the  Journal  for  color 
versions  of  Figs.  4  and  5]. 


(Garcia-Burillo  et  al.  1998)  show  evidence  of  an  inner 
triaxial  mass  distribution.  It  is  also  apparent  in  Figure  4  that 
the  spatial  isochrones  cover  only  a  fraction  of  the  active 
star-forming  regions  in  the  disks,  which  again  must  be 
caused  by  noncircular  kinematics.  However,  our  experimen¬ 
tation  with  rotation  curve  parameters  showed  that  to  first 
order,  our  simple  assumptions  for  the  rotation  curve  (flat, 
with  a  linear  inner  portion)  worked  best  in  matching  the  Ha 
morphology  of  the  arms. 

A  counterexample  is  shown  in  Figure  5  for  M51.  One 
previous  estimate  for  rco  in  this  galaxy  is  9.6  kpc  (Elmegreen 
et  al.  1992).  It  is  apparent  that  the  isochrones  do  not 
correspond  with  the  H  n  region  positions  as  cleanly  as  in 
Figure  4 a.  The  arms  show  a  smooth,  continuous  morphol¬ 
ogy  in  the  region  around  the  suggested  rco,  whereas  the 
predicted  isochrones  change  character  and  diverge. 


Table  2  shows  our  assumed  rotation  curve  parameters 
and  our  estimates  for  the  rco,  ILR,  and  OLR.  We  now 
compare  our  results  for  rco  with  those  in  the  literature,  in 
particular,  estimates  derived  from  kinematic  models.  There 
are  no  modeled  estimates  of  rco  for  NGC  695 1 . 

4.1.  M51 

Our  analysis  suggests  that  the  corotation  radius  is  twice 
as  large  as  most  previous  determinations  for  M51.  Both 
morphological  estimates  (e.g.,  Vogel  et  al.  1993;  Elmegreen 
et  al.  1992;  Knapen  et  al.  1992)  and  kinematic  model  esti¬ 
mates  (Garcia-Burrillo  et  al.  1993a)  found  rco  of  126"— 161", 
in  contrast  to  our  estimate  of  301".  We  note  that  M51  is  one 
of  the  few  examples  in  which  the  arm  definition  is  narrow 
enough  that  the  R-band  and  Ha  or  CO  observations  can 
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Fig.  4b 


clearly  distinguish  the  stellar  arms  from  the  new  star- 
formation  (Rand  &  Kulkarni  1993;  Garcia-Burrillo  et  al. 
1993b).  At  corotation,  the  youngest  regions  should  switch 
from  the  inside  edge  of  the  arm  to  the  outside;  yet  the  data 
shown  by  Rand  &  Kulkarni  (1990)  and  Garcia-Burrillo 
et  al.  (1993b)  show  no  evidence  of  this  transition  at  the 
smaller  rco  suggested  in  the  past. 


4.2.  M100 

For  Ml 00,  our  corotation  radius  is  larger  than  previous 
estimates,  but  by  a  smaller  margin.  Estimates  by  Garcia- 
Burrillo  et  al.  (1998)  and  Canzian  &  Allen  (1997),  based  on 
kinematic  modeling,  and  also  Elmegreen  et  al.  (1992)  all 
agree  fairly  well  with  the  range  of  82"— 1 14"  estimated  for  rco 
by  Sempere  et  al.  (1995).  Our  estimate  of  rco  ~  154"  is 
larger,  but  is  based  on  more  simplistic  morphological  and 
kinematic  assumptions;  most  of  the  above  studies  consider 
the  more  realistic,  bar-dominated  nature  of  the  rotation 


curve.  We  also  note  that  the  ends  of  the  arms  are  not  clearly 
delineated  in  this  galaxy. 

4.3.  NGC157 

Our  estimate  of  rco  ~  56"  agrees  well  with  estimates  in  the 
literature  based  on  detailed  kinematic  models.  Sempere  & 
Rozas  (1997)  find  rco  ~  50",  and  Fridman  et  al.  (2001)  find 
36"  <  rco  <  63"  by  two  different  methods.  Elmegreen  et  al. 
(1992)  also  estimate  corotation  to  be  at  56"  based  on 
morphological  arguments.  Spatial  isochrones  thus  appear 
to  offer  a  new  technique  for  constraining  the  corotation 
radius  for  a  first-order  approximation.  We  can  set  a  spatial 
correspondence  between  (1)  corotation  and  the  end  of  the 
nebular  spiral  arms,  (2)  isochrone  bunching  and  H  n  region 
clumps,  and  possibly  (3)  OLR  and  the  edge  of  the  star¬ 
forming  disk.  These  physical  coincidences  suggest  that  the 
location  of  corotation  has  been  identified  to  first  order. 
Thus,  while  some  authors  often  identified  the  ends  of  the 
star-forming  arms  with  the  corotation  resonance,  and 
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Fig.  4c 


others  with  the  OLR,  our  analysis  supports  the  former.  We 
also  note  that  in  wide-format,  blue  photographic  images 
available  via  the  NASA/IPAC  Extragalactic  Database 
(e.g.,  from  the  Sandage  &  Bedke  1988  atlas ),  there  is  no  evi¬ 
dence  in  any  of  our  galaxies  that  the  dust  lane  crosses  over 
from  the  inside  to  the  outside  of  the  arms  at  any  point  inside 
our  adopted  rco. 

5.  CONCLUSION 

To  better  understand  the  behavior  of  the  H  n  LF  and  its 
ability  to  probe  global  star  formation  properties,  we  have 
reexamined  its  relation  to  the  H  ii  region  size  distribution 
and  investigated  its  evolution  with  respect  to  density-wave 
triggering  in  spiral  galaxies. 


We  have  analytically  reevaluated  the  relationship 
between  the  H  n  LF  and  the  H  n  region  size  distribution. 
We  suggest  that  the  size  distribution  is  intrinsically  described 
as  a  power  law,  as  physically  related  to  the  H  n  LF,  rather 
than  its  conventional  description  by  an  exponential  relation. 
We  find  that  the  physical  relation  between  nebular  luminos¬ 
ity  and  size  implies  that  the  slopes  of  the  size  distribution  b 
and  the  luminosity  function  a  should  be  related  as 
b  =  2  —  3a.  The  normally  observed  value  a  =  2  thus  yields 
b  =  4,  which  we  find  to  be  consistent  with  observations.  We 
find  that  the  slope  transition  radius  log  RSSit  pc  ~  2. 1  pc  cor¬ 
responding  to  the  transition  luminosity  logLsat  ergs  s-1  ~ 
38.5  is  also  consistent  with  observations. 

We  have  developed  a  code  to  compute  spatial  iso¬ 
chrones  from  the  rotation  curve,  given  the  current  locus 
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Fig.  4  d 


of  the  spiral  density  waves.  Our  code  allows  quantitative 
comparison  of  the  nebular  spatial  distribution  with  the 
spatial  isochrones.  We  find  that  the  best  correspondence 
for  the  arm  H  n  regions  occurs  when  the  corotation 
radius  coincides  with  the  ends  of  the  star-forming  arms. 
This  value  for  rco  also  yields  a  coincidence  between  the 
inner  Lindblad  resonance  with  the  radial  onset  of  star 
formation  in  the  arms  for  three  of  the  four  galaxies. 
There  may  also  be  a  hint  that  the  outer  Lindblad  reso¬ 
nance  coincides  with  the  outer  edge  of  the  star-forming 
disk  in  at  least  one  galaxy.  The  spatial  isochrones  there¬ 


fore  appear  to  offer  a  new  technique  for  determining  the 
corotation  radius  and  principal  resonances. 

We  had  hoped  that  these  isochrones  would  isolate  subsets 
of  H  ii  regions  that  could  then  be  used  to  study  evolutionary 
effects  of  the  nebular  population.  However,  for  three  of  the 
four  galaxies,  the  isochrones  do  not  yield  adequate  spatial 
separation  over  the  nebular  evolution  timescale,  relative  to 
the  spatial  scatter  at  any  given  age.  Encouragingly,  we 
do  find  that  the  spatial  concentration  of  the  objects 
along  the  isochrones  is  consistent  with  the  overall  nebular 
evolutionary  timescale  in  these  galaxies.  However,  the  data 
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Fig.  5. — Same  as  Fig.  4a,  but  with  the  M51  corotation  at  9.6  kpc  suggested  by  Elmegreen  et  al.  (1992) 


do  suggest  periods  of  extended  star  formation  at  any  given 
location  that  blur  the  diagnostics.  Combined  with  uncer¬ 
tainties  in  the  kinematic  parameters,  this  greatly  limits  the 
ability  to  obtain  clear  quantitative  constraints  on  H  n  LF 
evolutionary  parameters  as  related  to  the  density-wave  pas¬ 
sage.  The  observed  rate  of  evolution,  however,  does  set 
limits  on  the  importance  of  this  process  in  producing  the 
observed  arm-interarm  variation  in  the  H  n  LF. 
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